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Abstract 

Three simple 7—, (7 + 3) — , 10— parametric trial functions for the H^^ molecular ion are presented. 
Each of them provides subsequently the most accurate approximation for the Born-Oppenheimer 
ground state energy among several-parametric trial functions. These trial functions are chosen 
following a criterion of physical adequacy and includes the electronic correlation in the exponential 
form ~ exp (7?"i2), where 7 is a variational parameter. The Born-Oppenheimer energy is found to 
be £' = —1.34034, —1.34073, —1.34159 a.u., respectively, for optimal equilateral triangular configu- 
ration of protons with the equilibrium interproton distance R = 1.65 a.u. The variational energy 
agrees in three significant digits (s.d.) with most accurate results available at present as well as 
for major expectation values. 
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I. INTRODUCTION 



The two-electron hydrogenic molecular ion Hjj" is one among the most abundant chemical 
compounds in the Universe. Its existence is of fundamental importance in chemistry and 
physics, in particular, because of its stability towards decay to H2 + p, the Hj]' ion is also 
a major proton donor in chemical reactions in interstellar space. The was discovered 
experimentally by J.J. Thomson in 1912 [1]. The system was very difficult for theoretical 
studies. Many theoretical methods were developed to study low-lying quantum states of 
this system. In particular, it became clear very quickly that interelectron correlation is 
of great importance and it should be included to the variational trial function explicitly 
which assure a faster convergence. This conclusion was similar to one drawn by James and 
Coolidge for the H2 molecule. Usually, the interelectron correlation was written in the form 
r^2 (Hylleraas [2| - James- Coolidge [3| form) or exp(— ar^g) (Gaussian form, see e.g. Ref. fl). 

Recently, Korobov \^ showed in explicit way that for the case of Helium atom the use of 
exponential form exp(— 7ri2) dramatically improves convergence and leads, in fact, to the 
most accurate results for the ground state energy for the Helium atom at present. Later on, 
it was shown that the similar use of exponential correlation exp(— 7ri2) for the H2 molecule 
allows to construct the most accurate trial function among few-parametric trial functions 
A hint why namely this ri2-dependence leads to the fast convergent results was given in 
In year 2006 an overwhelming discussion meeting took place in London, UK where different 
properties of the Hjj" ion and, in particular, various theoretical approaches to study the Hjj" 
ion were exposed (see [3]). 

The goal of this Note is to propose a simple, compact, easy-to-handle trial function 
depending exponentially on ri2 with few nonlinear parameters which leads to highly accurate 
Born-Oppenheimer ground state energy and major expectation values. We are not aware 
about previous studies of the H^ ion with trial functions involving ri2 in exponential form 
with a single exception js] where the Hj^ in linear configuration was explored. 

In this paper atomic units {h = e = rrif, = 1) are used throughout, albeit energies are 
given sometimes in Rydbergs. 



2 



II. THE H+ ION IN THE BORN-OPPENHEIMER APPROXIMATION 



The Hamiltonian which describes the ion under the assumption that the protons 
are infinitely massive (the Born-Oppenheimer approximation of zero order) and located at 
the vertices of an equilateral triangle of side R (see Fig. [T] for the geometrical setting and 
notations), is written as follows: 



K=A,B,C 



where pj = —iVj is the 3-vector of the momentum of the jth electron, the index k runs 
over protons A, B and C, rj^^ is the distance between the jth electron and the Kth proton, 
is the interelectron distance, and R is the interproton distance. 

It is a well established fact that the ground state of the Hjj' molecular ion is l^A'^, an elec- 
tronic spin-singlet state, with the three protons forming an equilateral triangle in the totally 
symmetric representation A\ of a Z^sh point symmetry [9]. Thus, the ground state elec- 
tronic wavefunction should be symmetric under permutations of the three indistinguishable 
protons. This ground state is the major focus of the present study. 

It is worth mentioning that the best theoretical value at the moment for the Born- 
Oppenheimer ground state energy is E = —1.34383562502 a.u. jlOl obtained with a basis 
of 1000 explicitly correlated spherical Gaussian functions with shifted centers. This value 
surpasses the previous record E = —1.343835624 a.u. by Cencek et al. which was obtained 
by using explicitly correlated Gaussian functions [ll[. 

III. VARIATIONAL METHOD 

The variational procedure is used as a method to explore the problem. The recipe of 
choosing the trial function is based on arguments of physical relevance, e.g. the trial function 



should support the symmetries of the system, has to reproduce the Cou 
and the asymptotic behavior at large distances adequately (see, e.g. [12 



omb singularities 



14|). In practice. 



the use of such trial functions implies the convergence of a special form of the perturbation 
theory where the variational energy is the sum of the first two terms. Let us remind the 
essentials of this perturbation theory (for details, see 12l-14|). Let us assume that our 



original Hamiltonian has a form % = —A + V. As a first step we choose a trial function 




FIG. 1: Geometrical setting for the hydrogen molecular ion Hg in equilateral triangular configu- 
ration. The three protons are located on the x-y plane forming an equilateral triangle with the 
origin of coordinates located at the geometrical center (circumcenter) of the triangle. 



^(triai) (^normalized to one) and find a potential for which such a trial function is an 

exact eigenfunction, i.e. Vtriai = A'tjj^^"'^^'^ / ip^^""''-\ with energy Etriai = 0. In a pure formal 
way we can construct a Hamiltonian Htriai = —A + Vtriai such that = 0. It can 

be easily shown that the variational energy 

is nothing but the first two terms in the perturbation theory where the unperturbed problem 
is given by l-Ltriai and the perturbation is the deviation of the original potential V from the 
trial potential Vtriau namely, Vperturhation = V — Vtriai- Eventually, we arrive at the formula 

Eyar Etriai ~l~ Ei(Vp^rturbation) i (2) 

here Ei{Vperturbation) = {i^^^""^^ \Vperturbation\'ip^^"''^^ ) IS the first energy correction in the per- 
turbation theory, where the unperturbed potential is Vtriai- It is worth noting that if the 
trial function is the Hartree-Fock function the resulting perturbation theory is nothing but 
the Moeller-Plesset perturbation theory (see, e.g. |l6i]. Section 15. IS)"*^. 

^ It is worth noting that the question about a convergence of the Moeller Plesset perturbation theory is not 
settled yet il7| 
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One of the criteria of convergence of the perturbation theory in Vperturbation = V — Vtriai 
is a requirement that the ratio \Vperturbation/V\ should not grow when r tends to infinity in 
any direction. If this ratio is bounded by a constant it should be less than one. In fact, it is 
a condition that the perturbation potential is subordinate with respect to the unperturbed 
potential. The value of this constant controls the rate of convergence - a smaller value of this 



constant leads to faster convergence [13|. Hence, the above condition gives a importance to 
the large-range behavior of the trial functions. In the physics language the above requirement 
means that the phenomenon of the Dyson's instability should not occur (for a discussion see 



12|) 



IV. CORRELATED TRIAL FUNCTION 



Among different forms to include explicit electronic correlation in the trial wave function 



for two-electron problems we mention three major approaches (see e.g. [18|): the linear in ri2, 
the gaussian exp(— ar^g) ^"^^ exponential exp(7ri2) terms. Among them, the factor exp(7ri2) 
only fulfills the adequacy requirements for a trial function described above. Thus, following 
the guidelines of Section III and the requirement of the convergence of the perturbation 
theory, we choose the trial function for the ground state in the following form: 

_ _|_ ^ ^ g-air'lA-"2»'lB-"3»"lC-"4''2A-"5»"2S-"6»'2C+7»'12 

pcrm{A, _B,C} 

where the sum runs over the permutations of the identical protons A, 5, C {S^ symmetry), 
and Pi2 is the operator which interchanges electrons (1 ^2). The variational parameters 
consist of non-linear parameters ai_6 and 7 which characterize the ( ant i) screening of the 
Coulomb charges. The interproton distance i?, see Fig.l is kept fixed. It is chosen to be 
equal R = R^q = 1.65a.u [10-11]. The function ([3]) is a symmetrized product of Is Slater 
type orbitals multiplied by the exponential correlation factor e'^^'^^. 

Calculations of the variational energy were performed using the minimization package MI- 
NUIT from CERN-LIB. Six-dimensional integrals which appear in the functional of energy 
were calculated numerically using a "state-of-the-art" dynamical partitioning procedure: the 



^ It is worth noting that this procedure for a selection of the trial function was applied successfully to a 
study of one-two-electron molecular systems in a magnetic field leading to the highly accurate results. 



Many of these results are the most accurate at the moment (see IJ] and 15| ) 
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domain of integration was divided into 972 subdomains following the profile of the integrand, 
in particular, separating out the domains with sharp changes of the integrand. Then each 
subdomain was integrated separately in parallel manner with controlled accuracy (for de- 



tails, see e.g. A realization of the routine required a lot of attention and care. During 

the minimization process the partitioning was permanently controlled and adjusted. Nu- 
merical integration of every subdomain was done with a relative accuracy of ~ 10"'^ — 10^^ 



depending on its complexity and re 



ative contribution using an adaptive routine based on 



an algorithm by Genz and Malik 19(| from R. Schiirer's HIntLib C++ multidimensional in- 



tegration library ( http : / / mint . sbg . ac . at /rudi/ J . Parallelization was reached using the MPI 



standard library MPICH. Computations were performed on a Linux cluster with 48 Xeon 
processors of 2.67 GHz each, and 12Gb total RAM plus extra processor serving as the master 
node. Total minimization process took about 1000 hours of wall clock time when a single 
call took about three minutes. For optimal values of parameters it took about 20 minutes 
(wall clock time) to compute a variational energy with relative accuracy 10~^. 

V. RESULTS 

In Table [T] we present the results for the ground state energy at interproton equilibrium 
distance of the molecular ion obtained by different researchers using different methods. 
In a clear way it is seen that the Born-Oppenheimer ground state energy obtained using the 
trial function ([3]) is the most accurate (the lowest) energy obtained with a few parametric 
functions. In particular, the trial function (|3]) gives a lower energy than the energies obtained 



with the explicitly correlated functions based on both Gaussians in 20| and linear in ri2 



2l| . when a relatively small number of terms with non-linear parameters is involved. The 
trial function (131) is more accurate than almost all(!) traditional CI calculations which were 

n 

performed before 1971 (see [221]) even including one of the largest set of 100 configurations 



24j |. In those CI calculations no explicit correlation was included. The variational energy 



obtained with ([3]) is even of comparable accuracy to the large CI calculations |23|, |2J] ^. 
Table [TTl shows the optimal values of the variational parameters in ([3]) . 

The list of major expectation values obtained with the trial function ([3]) and its com- 



^ for a list of 42 calculations of the ground state energy of H;^ in the period 1938-1992 see Ref.f22|, for a 



list of selected ab-initio calculations till 1995, see j^] 
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parison with results of other calculations is given in Table IIIII A reasonable agreement for 
expectation values is observed. In particular, for the expectation values (l/ri^), (a;^), {z^) 
and (r^) an agreement within ~ 1% with ours and all other calculations is observed, includ- 
ing ones obtained in the large C1SD-R12 calculations 2l|. Also, for the expectation value 
of (I/T12) we have an agreement in the first significant digit with other calculations being 
in closer agreement to the value obtained with the correlated Gaussian (unrestricted) wave- 
function with 15 terms [2^, while for the expectation values (l/ri^) and {z^) we observe 
an agreement with other calculations in 3 and 2 significant digits, respectively. These facts 
seem to indicate that the presented expectation values are very accurate, corroborating the 
quality of the trial function (|2]) giving 2-3 s.d. correctly. Perhaps, it is worth noting that in 
absence of any criteria about accuracy of the obtained expectation values we can only note 
about agreement of them obtained in different approaches. 

VI. CONCLUSION 

We presented a simple and compact 7-parametric variational trial function together 
with its possible natural generalization by addition of the Heitler-London (HL) type func- 
tion. This function already provides surprisingly accurate Born-Oppenheimer energy for the 
ground state of such a complicated molecular system Hjj'. It is chosen following a criterion of 
physical adequacy which suggests to take the electronic correlation in the form ~ exp (7ri2) 
where 7 is a variational parameter. The minimum energy is found iohe E = — 1.34034 a.u. 
at an equilibrium interproton distance R = 1.65 a.u. This result for the energy is the most 
accurate among the values obtained with several parametric trial functions. In particular, 
it is more accurate than the energies obtained with the explicitly correlated approaches of 



Ref. 2l|] (linear in ri2) and that of Ref.[20| (Gaussian in ri2), when a relatively small number 
of terms and non-linear parameters are involved. 

n n 

In a spirit of the approach presented in [6| (see also [5[) the trial function ([3]) can be 
modified by adding similar function, in particular, of the Heitler-London type: 

— Q-a{riA-riB-ric-r2A-r2B-r2c)+^ri2 ^^■^ 

where a, 7 are parameters. The function alone gives a dominant contribution to small 
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E (a.u.) 


R (a.u.) 


method 


reference 


-1.339 7 


1.66 


CI-GTO, > 120 configs 


[25] (1970) 


-1.306 29 


1.65 


GG, 3 terms, 5 non-linear params 


[201 (1973) 


-1.327 25 


1.65 


GG, 6 terms, 7 non-linear params 




-1.331 47 


1.65 


GG, 10 terms, 9 non-linear params 




-1.332 29 


1.65 


GG, 15 terms, 11 non-linear params 




-1.334 382 


1.65 


R12, 10s basis set 


[21J (1993) 


-1.334 632 


1.65 


R12, 30s basis set 




-1.340 34 


1.65 


7-Parametric Trial Function ([3]) 


present 


-1.340 5 


1.6405 


CI -GTO, 48 configs 


[23] (1971) 


-1.340 5 


1.65 


CI -STO, 100 configs 


[24] (1971) 


-1.340 73 


1.65 


(7-|-3)-Parametric Trial Function ([5]) 


present 


-1.341 59 


1.65 


10-Parametric Trial Function ([5]) 


present 


-1.342 72 


1.65041 


CI-GTO, 108 terms 


[27] (1985) 


-1.342 784 


1.6504 


CI-GTO, 8s3pld/[6s3pld] basis set 


[28] (1978) 


-1.343 40 


1.6504 


CI-GTO, I0s4p2d basis set 


[29] (1990) 


-1.343 822 


1.65 


CI-GTO, 700 terms 


[30] (1990) 


-1.342 03 


1.6504 


CI with rl2, 36 configs 


[31] (1982) 


-1.343 422 


1.6504 


CI with rl2, 192 configs 


[32] (1984) 


-1.343 500 


1.6504 


CI with rl2, 13s3p/[10s2p] basis set 


[33] (1988) 


-1.343 828 


1.65 


CI with rl2, 13s5p3d basis set. 


[34] (1990) 


-1.343 835 


1.65 


R12, 30s20pl2d9f basis set 


[21] (1993) 


-1.343 35 


1.65 


GG, 15 terms, 135 non-linear params 


[20] (1973) 


-1.343 835 624 


1.65 


GG, 600 terms 


[IJJ (1995) 


-1.343 835 625 02 


1.65 


ECSG, 1000 terms 


[10] (2009) 



TABLE I: A selection of the calculations for the Born-Oppenheimer ground state energy at equi- 
librium distance of Hj^. Record calculations of Ref. |lo| (2009) and Ref. 11] (1995). CI denotes 



Configuration Interaction, STO - Slater Type Orbitals, GTO - Gaussian Type Orbitals, GG - 
correlated Gaussians (Gaussian Geminals), R12 - the CI calculation augmented by terms linear in 
ri2, ECSG - Explicitly Correlated Spherical Gaussian functions. Trial Function ([5]) with the 

parameters of ipo kept fixed and equal to ones found for ([SD, Trial Function ([5]) with all 10 

■ ■ , 8 

parameters optimized. 



E{Ry) 


«! a2 Q!3 04 05 ag 7 


A a J 


-2.680 7 


-0.00353 0.18548 1.4245 1.0471 0.15082 0.58912 0.21632 




-2.681 4 




-0.03000 0.47517 0.76398 


-2.683 2 


-0.00294 0.21022 1.3849 1.0199 0.17103 0.59084 0.26044 


-0.51154 0.59589 0.86229 



TABLE II: The ground state energy of at Reg = 1.65 a.u. and the non-hnear variational param- 
eters in [a.u.]^^ corresponding to the trial function ([3]), to the trial function ([5]) with parameters 
corresponding to ipo fixed and to ([5]) with 10 optimized parameters. 



interproton distances. Taking a linear superposition with 



(5) 



and making minimization with respect to parameters A, a and 7 only (see Table [ITl) give 
an essential improvement in the energy (see Table [T]). In particular, this function, which 
contains (7 + 3) variational parameters, allows us to get more accurate result for energy 



than one obtained in [2^ within CI-STO with 100 configurations. 

Releasing all 10 parameters in ([5]) (see Table [Tll) we obtain further improved result (see 
Tabled]), although we are still unable to reproduce the fourth significant figure in the energy. 
However, the obtained energy is among the thirteen the most accurate variational results 
ever calculated so far. It is slightly worse then one obtained in j3l| based on CI with rl2 
method with 36 configurations. The expectation values in Table IIIII gradually change with 
move from one Ansatz to another seemingly demonstrating a convergence. It seems evident 
that taking a linear superposition of two (or more) functions ([3]) instead of ([5]) will improve 
essentially the variational energies. It will be done elsewhere. Undoubtfully, trial functions 
([3]), can be used to study potential energy surface. It is worth emphasizing that the 
main attraction of functions ([3]), ([S]) is their compactness. 

The function ([3]) can be easily modified for a study of spin-triplet states and as well 
as the low-lying states with non-vanishing magnetic quantum number. A generalization to 
more-than-two electron molecular systems seems also straightforward. 
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Expectation Value 


iriai runction / (|oj)^ / (O*- ' 


Others 




2.0032 






z.OOio 






1 nnQ 1 

i.yyoi 




(IA12) 


0.6315 


0.59415" 




(J.ooU4 


O.DZDOD 




O.DoOz 




(IAia) 


0.8548 


0.85519'= 




0.8548 


0.8553 




0.8549 




lJ2\ l„ 2\ 


0.77ii 


0.75818 




0.7703 


0.7595'' 




7fififi 


7591 3^= 






0.75968 






0.7605 


/^2\ 


0.5399 


0.54802 




0.5367 


0.5451'' 




0.5337 


54085*= 






0.54179 






0.5396^= 


(r ) 


z.08zz 


2.06442 




2.0773 


2.0640'' 




2.0669 


2.05911"= 






2.06114"^ 



TABLE III: Expectation values (in a.u.) for the Hg" ion in its ground state obtained with the 
trial functions ([3]) and (jS])*^*'**). Corresponding results obtained with other methods are displayed 
for comparison. Coordinates x, y, z and r are measured from the center of the equilateral triangle 

m 



formed by protons. Ref. [23l] CI-43; '' CI wavefuncion (I) in Ref. 
(unrestricted) wavefunction with 15 terms in Ref. 2C|; CI wavefunction in Ref. 



wavefunction with the 10s8p6d4/ basis set in Ref. 



Correlated Gaussian 
261; CISD-R12 
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